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Executive  Summary 

All  industrial  and  military  organizations  rely  on  the  effective  use  of  repairable  systems.  In  many  cases, 
these  organizations  rely  on  sets  of  identical  systems  that  are  required  to  perform  sets  of  missions.  Often, 
sufficient  maintenance  resources  aren’t  available  to  keep  all  systems  within  a  set  in  top  condition  at  all 
times;  therefore,  a  maintenance  manager  must  decide  how  to  allocate  the  available  resources.  This 
allocation  falls  within  the  domain  of  selective  maintenance.  Selective  maintenance  is  defined  as  the 
process  of  identifying  the  subset  of  maintenance  activities  to  perform  from  a  set  of  desired  maintenance 
actions.  Although  the  modeling  of  repairable  equipment  has  been  studied  extensively,  traditional  studies 
tend  to  focus  on  a  single  system  and  ignore  the  mission  profile  of  the  system.  For  the  United  States  Air 
Force  (USAF),  these  limitations  prevent  current  models  from  providing  meaningful  guidance  relative  to 
maintenance  planning  and  sortie  scheduling. 

The  objective  of  this  project  is  to  investigate  the  use  of  a  mathematical  modeling  methodology  for 
integrating  maintenance  planning  and  sortie  scheduling  issues.  First,  the  relevant  research  literature  for 
both  selective  maintenance  and  fleet  assignment  is  presented.  The  selective  maintenance  literature  is 
limited  in  that  current  models  only  consider  a  single  system,  and  most  of  the  fleet  assignment  literature  is 
motivated  by  the  commercial  airline  industry,  which  focuses  on  meeting  Federal  Aviation  Administration 
(FAA)  maintenance  requirements. 

Next,  background  research  is  presented,  which  extends  a  current  selective  maintenance  model  to 
incorporate  sets  of  systems.  Here  a  selective  maintenance  model  for  a  set  of  systems  that  must  execute  a 
set  of  missions  with  system  maintenance  performed  only  between  missions  is  defined.  Originally,  the 
scenario  evaluated  in  the  background  research  was  done  using  a  total  enumeration  strategy.  Although  the 
enumeration  strategy  guarantees  an  optimal  solution,  scenarios  involving  more  than  four  systems  result  in 
lengthy  computational  times;  therefore,  a  genetic  algorithm  (GA)  was  developed  to  solve  larger  scenarios. 
This  background  research  and  the  development  of  the  GA  contains  key  concepts  used  in  the  formulation 
and  solution  procedures  of  both  the  static  and  dynamic  optimization  models  that  integrate  aircraft 
assignment  and  selective  maintenance. 

The  model  developed  in  the  background  research  does  not  address  mission  assignment,  because  all 
missions  in  the  upcoming  set  are  identical,  hi  an  extension  to  this  model,  we  formulate  a  mathematical 
optimization  model  that  integrates  system  assignment  (given  a  relatively  static  mission  profile)  and 
selective  maintenance  decision-making.  Specifically,  we  consider  upcoming  missions  that  are  not 
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identical  and  define  a  parameter  to  quantify  the  “difficulty”  of  each  upcoming  mission.  We  develop  a 
solution  procedure  for  performing  this  optimization  and  study  the  behavior  of  the  model  using  numerical 
examples.  Because  of  the  complexity  of  the  optimization  model,  we  utilize  a  GA  to  perform  the  required 
optimization. 

Finally,  we  formulate  a  more  complex  optimization  model  that  addresses  a  more  dynamic  mission  profile. 
Specifically,  missions  start  and  end  at  different  times,  and  maintenance  and  scheduling  decisions  are 
made  over  a  series  of  time  “buckets.”  We  consider  a  planning  horizon  such  that  each  system  in  the  set 
returns  from  its  previous  mission  and  begins  its  future  mission;  however,  no  system  returns  from  its  future 
before  the  end  of  the  planning  horizon.  We  again  utilize  GAs  to  perform  the  required  optimization. 
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1  Introduction 


All  military  organizations  depend  upon  the  reliable  performance  of  repairable  systems  for  the  successful 
completion  of  missions.  Although  the  use  of  mathematical  modeling  for  the  purpose  of  modeling 
repairable  systems  and  designing  optimal  maintenance  policies  for  these  systems  has  received  an 
extensive  amount  of  attention  in  the  literature,  this  research  is  limited  in  two  key  ways.  First,  studies  in 
the  literature  tend  to  focus  on  a  single  system  rather  than  a  set  of  systems  used  by  the  organization. 
Second,  these  studies  ignore  the  mission  profile  of  the  system,  which  prevents  the  modeler  from 
considering  important  maintenance  strategies  including: 

♦  Performing  maintenance  during  scheduled  downtime 

♦  Delaying  maintenance  to  execute  a  critical  mission 

For  the  USAF,  these  limitations  are  too  severe  to  provide  meaningful  guidance  relative  to  fleet 
maintenance  planning.  For  a  single  aircraft,  maintenance  actions  should  be  planned  relative  to  its  mission 
schedule,  and  sortie  scheduling  decisions  should  be  managed  with  considerations  for  aircraft 
maintenance.  Given  that  aircraft  within  the  USAF  share  maintenance  resources  (spares,  labor,  etc.)  and 
perform  missions  in  groups,  the  integration  of  sortie  scheduling  and  maintenance  planning  can  become 
quite  complex. 

The  objective  of  this  project  is  to  investigate  the  use  of  a  mathematical  modeling  methodology  for 
integrating  maintenance  planning  and  sortie  scheduling  issues.  Achieving  the  objective  of  the  project 
requires  the  completion  of  five  key  activities. 

First,  we  define  a  hypothetical  set  of  systems  that  could  correspond  to  aircraft,  subsystems  within  an 
aircraft,  or  subsystems  within  subsystems.  This  definition  includes  the  number  of  systems,  the  mission 
profile  and  the  constrained  maintenance  resources. 

Second,  we  formulate  a  mathematical  optimization  model  that  integrates  system  assignment  (given  a 
relatively  static  mission  profile)  and  selective  maintenance  decision-making.  Note  that  selective 
maintenance  refers  to  the  process  of  identifying  the  subset  of  actions  to  perform  from  a  set  of  desirable 
maintenance  actions. 
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Third,  we  develop  a  solution  procedure  for  performing  this  optimization  and  study  the  behavior  of  the 
model  using  numerical  examples.  Because  of  the  complexity  of  the  optimization  model,  we  use  a  GA  to 
perform  the  required  optimization. 

Fourth,  we  formulate  a  more  complex  optimization  model  that  addresses  a  more  dynamic  mission  profile. 
And  finally,  we  again  utilize  GAs  to  perform  the  required  optimization. 

The  remainder  of  this  report  is  organized  as  follows.  Section  2  summarizes  the  relevant  research 
literature.  Section  3  details  a  description  of  background  research  critical  to  this  project.  Section  4 
formulates  the  optimization  model  for  the  static  scenario  and  describes  the  GA-based  solution  approach. 
Section  5  contains  a  similar  presentation  for  the  dynamic  scenario.  Finally,  Section  6  describes  future 
opportunities. 
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2  Research  Literature  Review 


This  project  builds  upon  the  body  of  knowledge  in  both  selective  maintenance  and  aircraft  fleet 
assignment;  therefore,  in  this  section,  we  summarize  the  relevant  literature  in  both  areas. 

2.1  Selective  Maintenance  Research 

Selective  maintenance  falls  under  the  domain  of  maintenance  modeling  and  optimization.  Fortunately,  the 
use  of  mathematical  modeling  for  the  purpose  of  modeling  repairable  systems  and  designing  optimal 
maintenance  policies  for  these  systems  has  received  an  extensive  amount  of  attention  in  the  literature  [8, 
10,  14,  15,  18,  19,  20].  The  original  study  in  selective  maintenance  was  performed  by  Rice  et  al.  [16], 
They  define  a  system  that  must  complete  a  series  of  missions  where  maintenance  is  performed  only 
during  finite  breaks  between  missions.  Due  to  the  limited  maintenance  resource  (in  this  case,  time),  it  may 
not  be  possible  to  repair  all  failed  components  before  the  next  mission.  A  nonlinear,  discrete  selective 
maintenance  optimization  model  is  developed  that  is  designed  to  maximize  system  reliability  for  the  next 
mission.  The  numbers  of  components  to  be  repaired  are  the  decision  variables,  and  the  limitation  on 
maintenance  time  serves  as  the  primary  functional  constraint.  Due  to  the  complexity  of  the  model,  total 
enumeration  is  the  recommended  solution  procedure.  Given  that  total  enumeration  is  ineffective  for  large 
scenarios,  a  heuristic  selective  maintenance  procedure  is  developed. 

Cassady  et  al.  [3,  4]  extend  the  work  of  Rice  et  al.  [16]  in  several  ways.  First,  a  more  complex  system  is 
defined.  Specifically,  a  system  is  comprised  of  independent  subsystems  connected  in  series  with  the 
individual  components  in  each  subsystem  connected  in  some  fashion.  Next,  the  selective  maintenance 
model  is  extended  to  consider  the  case  where  both  time  and  cost  are  constrained.  This  leads  to  the 
development  of  three  different  selective  maintenance  models.  These  models  include  maximizing  system 
reliability  subject  to  both  time  and  cost  constraints;  minimizing  system  repair  costs  subject  to  a  time 
constraint  and  a  minimum  required  reliability  level;  and  minimizing  total  repair  time  subject  to  both  cost 
and  reliability  constraints. 

Cassady  et  al  [5]  build  upon  the  work  of  Rice  et  al.  [16]  in  two  other  ways.  First,  system  components  are 
assumed  to  have  Weibull  life  distributions.  This  assumption  permits  systems  to  experience  an  increasing 
failure  rate  (IFR)  and  requires  monitoring  of  the  age  of  components.  Second,  the  selective  maintenance 
model  is  formulated  to  include  three  maintenance  actions:  minimal  repair  of  failed  components, 
replacement  of  failed  components,  and  replacement  of  functioning  components  (preventive  maintenance). 
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Chen  et  al.  [7]  expand  the  work  of  Rice  et  al.  [16]  and  Cassady  et  al.  [3]  by  considering  systems  in  which 
each  component  and  the  system  may  be  in  K  +  l  possible  states,  0,  1,  ...  ,  K.  They  use  an  optimization 
model  to  minimize  the  total  cost  of  maintenance  activities  subject  to  a  minimum  required  system 
reliability. 

2.2  Fleet  Assignment  Research 

Several  studies  in  the  literature  address  aircraft  fleet  assignment.  Abara  [1]  formulates  a  mixed-integer 
programming  fleet  assignment  model.  The  objective  function  of  the  model  may  maximize  profit, 
minimize  cost,  or  optimize  the  use  of  a  specific  type  of  aircraft.  He  assumes  a  certain  amount  of  “stand¬ 
by”  time  for  routine  maintenance  between  arrivals  and  departures.  Due  to  the  computational  difficulties 
associated  with  modeling  fleet  assignment  as  a  mixed-integer  program,  Hane  et  al.  [12]  investigate  the 
effect  of  several  optimization  improvement  methodologies  on  computational  time.  Rushmeier  et  al.  [17] 
improved  upon  the  work  of  Abara  [1]  and  Hane  et  al.  [12]  by  adding  “sit  activity  arcs,”  w'hich  can  be  used 
to  distinguish  between  aircraft  on  the  ground.  This  extension  could  allow  for  a  more  meaningful  analysis 
on  maintenance  and  crew  resource  utilizations. 

Lam  [13]  provides  a  brief  history  of  commercial  aircraft  maintenance.  He  highlights  the  effects  of 
economics,  technology,  and  regulation  of  the  FAA  on  aircraft  maintenance.  The  different  types  of  aircraft 
maintenance  required  by  the  FAA  are  defined  and  compared  by  the  amount  of  labor  hours,  maintenance 
hours,  and  the  maintenance  interval  hours  required  by  an  average  aircraft.  The  actual  maintenance  tasks 
performed  on  an  aircraft  are  also  defined  for  each  type  of  maintenance  check.  Barnhart  et  al.  [2]  develop  a 
modeling  methodology  for  integrating  fleet  assignment  and  aircraft  routing.  They  develop  an  example 
using  data  from  commercial  airlines  and  determine  that  the  proposed  model  successfully  addresses  both 
fleet  assignment  and  aircraft  routing.  Feo  et  al.  [11]  develop  a  model  that  can  be  used  to  locate 
maintenance  stations  and  to  develop  flight  schedules  that  meet  maintenance  requirements.  Specifically, 
the  model  seeks  to  optimize  the  number  of  required  maintenance  facilities  to  meet  FAA  maintenance 
regulations.  The  network  model  is  solved  using  a  heuristic  and  demonstrated  using  data  obtained  from 
American  Airlines.  The  model  is  applied  to  the  American  Airlines  Boeing  727  fleet,  and  the  results  show 
significant  improvement  over  other  aircraft  scheduling  models.  Similar  results  are  obtained  when  the 
model  is  applied  to  the  Super  80  and  DC-10  fleets.  Clarke  et  al.  [9]  expand  on  the  model  developed  by 
Hane  et  al.  [12]  by  adding  maintenance-  and  crew'-related  constraints  to  the  fleet  assignment  model.  The 
maintenance  constraints  address  both  short-  and  long-duration  maintenance. 
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3  Background  Research 

Previously  published  research  in  the  area  of  selective  maintenance  focuses  on  a  single  system.  In  an 
unpublished  research  effort  preceding  this  project,  Cassady  and  Schneider  [6]  formulate  an  optimization 
model  to  extend  the  work  of  Rice  et  al.  [16]  by  defining  a  selective  maintenance  model  for  a  set  of 
systems  that  must  perform  a  set  of  missions  with  system  maintenance  performed  only  between  missions. 
This  background  research  contains  key  concepts  used  in  the  formulation  of  both  the  static  and  dynamic 
optimization  models  that  integrate  aircraft  assignment  and  selective  maintenance. 

3.1  Hypothetical  Set  of  Systems 

Consider  a  set  of  q  independent  and  identical  systems.  Each  system  is  comprised  of  m  independent 
subsystems  connected  in  series,  and  each  subsystem,  say  subsystem  j,  contains  n,  independent  and 
identical  copies  of  a  constant  failure  rate  (CFR)  component  connected  in  parallel.  A  graphical 
representation  (reliability  block  diagram)  of  an  example  set  of  systems  is  shown  in  Figure  1.  Note  that 
each  component  in  Figure  1  is  labeled  “/,  j,  k”  where  i  denotes  the  system  number,  j  denotes  the 
subsystem  number  and  k  denotes  the  component  number. 

The  set  of  systems  is  required  to  perform  sequential  sets  of  missions  that  originate  from  and  return  to  a 
common  base  of  operation.  Each  set  of  missions  contains  q  independent  and  identical  missions  with 
common  start  times  and  durations.  At  any  time,  a  specific  component,  subsystem  or  system  is  in  one  of 
two  states:  functioning  (1)  or  failed  (0).  Components  only  fail  during  missions,  and  failed  components  can 
be  repaired  only  during  the  breaks  between  sets  of  missions.  Note  that  all  maintenance  is  performed  at  the 
base. 
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Figure  1 :  Example  Set  of  Systems 

3.2  System  Performance 

Suppose  that  a  set  of  missions  has  just  ended  and  all  systems  have  returned  to  their  base  of  operation  and 
maintenance.  Note  that  component  ijk  refers  to  component  k  in  subsystem  j  of  system  /,  and  note  that 
subsystem  ij  refers  to  subsystem  j  in  system  i.  Let  denote  the  status  of  component  ijk,  let  Yy  denote  the 
status  of  subsystem  ij\  and  let  Y{  denote  the  status  of  system  i.  Since  each  subsystem  is  a  parallel 
arrangement  of  its  components,  then 

Equation  3.1 

Since  each  system  is  a  series  arrangement  of  its  subsystems,  then 

m 

?-nr. 

.M 

Equation  3.2 

f  ■ 
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Note  that  if  Yy  =  0  for  some  j  =  1 ,  2, . . .  ,  m,  then  Y,  =  0  and  system  i  did  not  successfully  complete  the  last 
mission. 

Prior  to  the  next  set  of  missions,  maintenance  (repair)  can  be  performed  on  some  or  all  of  the  failed 
components.  Let  Xyk  denote  the  status  of  component  ijk  at  the  start  of  the  next  mission,  let  Xy  denote  the 
status  of  subsystem  ij  at  the  start  of  the  next  mission,  and  let  Xj  denote  the  status  of  system  /  at  the  start  of 
the  next  mission.  Note  that 

k  =  1 

Equation  3.3 


and 


m 

>1 

Equation  3.4 

Furthermore,  note  that  if  Yjjk  =  0  and  XJk  =  1 ,  then  component  ijk  was  repaired  during  the  break  between 
missions.  Also,  note  that  if Xg  =  0  for  some/  =  1,2,  ...  ,  nu  then X,  =  0  and  system  /  is  unable  to  undertake 
its  next  mission. 

Clearly,  the  Xuk  values  have  a  direct  effect  on  the  status  of  components,  subsystems  and  systems  at  the 
end  of  the  upcoming  mission.  Let  Yl  denote  the  status  of  component  ijk  at  the  end  of  the  next  mission, 

let  YJ  denote  the  status  of  subsystem  ij  at  the  end  of  the  next  mission,  and  let  Yj  denote  the  status  of 
system  i  at  the  end  of  the  next  mission.  Note  that 


and 


*=i 

Equation  3.5 
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m 

./'=! 

Equation  3.6 

Also,  note  that  if  Xyk  =  0,  then  Yjk  =  0 .  If  Xyu  —  1  then  YtJk  is  a  random  variable.  Assuming  Xt  =  1 ,  then 
both  Y~  for  ally  =  1,  2,  ...  ,  m  and  Y+  are  random  variables.  Furthermore,  note  that  if  Y*  =  0  for  some  j 
-  1, 2, ... ,  m,  then  Y *  =  0  and  system  /  did  not  successfully  complete  the  next  mission. 

One  key  factor  in  making  maintenance  decisions  is  the  probability  that  each  system  will  successfully 
complete  its  next  mission.  We  assume  that  all  components  have  constant  failure  rates,  that  is,  the 
probability  that  a  component  successfully  completes  a  mission  depends  only  upon  its  status  at  the 
beginning  of  the  mission  and  the  length  of  the  mission.  Since  all  missions  have  identical  length,  we  ignore 
mission  length  in  our  analysis.  Let 

0=Prfe=lk„=l) 


Equation  3.7 

Note  that  this  probability  is  only  indexed  on  j  because  systems  are  identical  and  the  components  within  a 
subsystem  are  identical.  The  reliability  of  component  ijk  is  defined  to  be  the  probability  that  component 
ijk  is  functioning  at  the  end  of  the  next  mission  and  given  by 

=  Prfe  =  i)=  ■ 

Equation  3.8 

Likewise,  the  reliability  of  subsystem  ij  is  defined  to  be  the  probability  that  subsystem  ij  is  functioning  at 
the  end  of  the  next  mission  and  given  by 

*=1  k= 1 

Equation  3.9 
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Since  the  components  within  a  subsystem  are  identical.  Equation  3.9  can  be  simplified.  Let  by  denote  the 
number  of  functioning  components  in  subsystem  ij  at  the  beginning  of  the  next  mission.  Note  that 

»,=£** 

k= 1 

Equation  3.10 

Then 

Equation  3.11 

Finally,  the  reliability  of  system  i  is  defined  to  be  the  probability  that  system  i  is  functioning  at  the  end  of 
the  next  mission  and  given  by 

m 

>i 

Equation  3.12 

3.3  Selective  Maintenance  Model  for  a  Set  of  Systems 

Recall  that  we  consider  the  case  in  which  a  set  of  missions  has  just  ended  and  all  systems  have  returned  to 
their  base  of  operation  and  maintenance.  The  base  possesses  the  technological  capability  to  repair  any 
failed  component  and  ideally,  all  failed  components  would  be  repaired  prior  to  the  beginning  of  the  next 
mission;  however,  maintenance  resource  limitations  may  prevent  the  repair  of  all  failed  components. 

Let  cijj  denote  the  number  of  failed  components  in  subsystem  ij,  and  note  that  ay  is  given  by 

“,=£(1-4) 

A=] 

Equation  3.13 

Furthermore,  the  total  number  of  failed  components  of  type  j  is  given  by 
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</ 

a-j='Zau 

M 

Equation  3.14 

Each  repair  consumes  a  fixed  amount  of  each  of  s  maintenance  resources  (for  example,  labor,  equipment, 
spare  parts).  Let  cty  denote  the  amount  of  resource  /  consumed  by  repairing  a  component  of  type  j,  and 
suppose  fa  denotes  the  amount  of  resource  /  available  during  a  single  break.  If 

m 

YuaHav  - A 
>1 

Equation  3.15 

for  all  /=  1,  2,  ...  ,  s,  then  all  failed  components  may  be  repaired  prior  to  the  next  mission.  Otherwise,  a 
method  is  needed  to  decide  which  failed  components  should  be  repaired  and  which  components  must 
remain  in  a  failed  condition.  This  decision-making  process  is  the  selective  maintenance  scenario  that  we 
consider. 

Let  djj  denote  the  number  of  failed  components  in  subsystem  ij  to  be  repaired  prior  to  the  beginning  of  the 
next  set  of  missions.  Then 

4  =£(*»» -4) 

k= 1 

Equation  3.16 

The  total  number  of  failed  components  of  type  j  repaired  prior  to  the  next  set  of  missions  is  given  by 

i=i 

Equation  3.17 

To  address  the  selective  maintenance  issue,  we  formulate  a  nonlinear,  discrete  optimization  model  where 
the  djj  values  serve  as  the  decision  variables. 
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The  first  constraint  on  these  decision  variables  is  that  they  be  integer-valued.  Second,  the  number  of 
repairs  is  limited  to  the  number  of  failed  components,  that  is, 

0  <  dy  <  a,y 

Equation  3.18 

Third,  the  number  of  repairs  is  limited  by  the  available  maintenance  resources,  that  is, 

m 

7=1 

Equation  3.19 

for  all  /  =  1 ,  2,  . . .  ,  s.  The  objective  in  choosing  values  for  the  decision  variables  is  to  maximize  the 
probability  that  all  upcoming  missions  are  completed  successfully.  This  probability  is  referred  to  as 
overall  reliability  and  given  by 


i~\  i-\  7=1 


where 


Equation  3.20 


bu  =  ni 


aii+di.i 


Equation  3.21 

The  resulting  optimization  model  is  given  by 
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Maximize 


P: 


9  »>  ,  , 

W  j-\ 


subject  to 


<7  m 

'YJEjajidij  ^  Pi 


/=)  7=1 


7=1,2,  ... 


A’ 


0  <  di}  <  a0 ,  integer  i=\,2,...,q 

j=],2,...,m 

Equation  3.22 

As  a  numerical  example,  consider  the  system  presented  in  Figure  1  having  q  =  2  and  m  =  3.  Two  (s  =  2) 
limited  maintenance  resources  limit  repair  activities  during  the  break  between  mission  sets.  Note  that  /3\  = 
16  and  =  10.  The  remaining  parameters  for  the  set  of  systems  are  given  in  Table  3.1.  Suppose  all 
systems  have  just  returned  from  a  set  of  missions.  The  status  of  the  set  of  systems  is  summarized  in  Table 
3.2  and  Figure  2.  Note  that  both  systems  successfully  completed  the  last  mission.  Repairing  all  the  failed 
components  would  require  18  units  of  resource  1  and  11  units  of  resource  2.  Therefore,  resource  1  is 
limited  and  we  can  use  the  selective  maintenance  model  P  to  determine  which  components  to  repair  to 
maximize  the  overall  reliability.  By  enumerating  all  feasible  solutions,  we  identify  the  optimal  selective 
maintenance  solution  for  this  scenario.  This  solution  is  summarized  in  Table  3.3,  and  the  revised  status  of 
the  set  of  systems  based  on  the  solution  is  shown  in  Figure  3.  If  no  maintenance  were  performed,  the 
overall  reliability  would  be  0.5688,  and  if  sufficient  maintenance  resources  were  available  to  repair  all 
failed  components,  the  overall  reliability  would  be  0.9665. 


Table  3.1:  Model  Parameters  for  Example  Set  of  Systems 


Subsystem 

D 

«yi 

Op 

1 

2 

0.90 

2 

01 

2 

3 

0.85 

3 

■ 

3 

2 

0.94 

2 

m 
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Table  3.2:  Failed  Components  for  Example  Set  of  Systems 


System  (q) 

Subsystem  (0 

1 

2 

3 

1 

an  =  1 

ai2  -  2 

ai3  =  1 

2 

a2i  =  0 

a22  “  2 

a23  =  1 

Figure  2:  Status  Prior  to  Mission 


Table  3.3:  Optimal  Selective  Maintenance  Actions 


System 

cfo 

m 

da 

Reliability 

Fleet  Reliability 

1 

1 

2 

■ 

0.9831 

0.9479 

2 

o 

1 

■ 

0.9642 
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Figure  3:  Status  after  Maintenance 

3 A  Solution  Procedures 

For  small  scenarios  such  as  the  one  presented,  Visual  Basic®  code  that  enumerates  all  feasible  solutions  is 
executed  in  a  macro  within  a  Microsoft®  Excel  spreadsheet.  The  code  generates  a  solution  (specifically, 
the  djj  values),  and  the  solution’s  feasibility  is  checked  against  the  model  constraints.  If  the  solution  is 
infeasible,  it  is  disregarded;  however,  if  a  feasible  solution  is  generated,  the  resulting  overall  reliability  is 
calculated.  If  the  overall  reliability  of  the  solution  is  equal  to  or  greater  than  the  largest  reliability 
computed  previously,  that  solution  (and  its  corresponding  overall  reliability)  is  written  to  a  file  and  a  new 
solution  is  generated.  After  generating  all  possible  solutions,  the  optimal  solution  (or  all  optimal  solutions 
in  the  case  of  a  tie)  is  output  to  an  Excel  worksheet.  This  enumeration  code  is  able  to  handle  any  scenario 
size  and  guarantees  the  identification  of  all  optimal  solutions.  Also,  solution  feasibility  checks  and 
computations  are  done  “on-the-fly,”  eliminating  the  excessive  use  of  computer  memory. 

Although  the  enumeration  code  guarantees  an  optimal  solution,  scenarios  involving  more  than  four 
systems  result  in  lengthy  computational  times.  Therefore,  the  use  of  a  search-based  heuristic  based  on  the 
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concepts  of  GAs  is  applied  to  larger  scenarios.  GAs  are  computationally  intense  search  methods 
originally  applied  in  the  area  of  artificial  intelligence.  Over  the  past  decade,  increases  in  computer  speeds 
have  made  this  directed  trial-and-error  method  increasingly  popular.  GAs  operate  on  the  principle  of 
“survival  of  the  fittest.”  In  general,  GA  keeps  track  of  a  fixed  population  of  candidate  solutions  called 
chromosomes.  Each  element,  or  gene,  within  a  chromosome  represents  some  portion  of  the  solution.  The 
solutions  are  then  ranked  using  a  fitness  function  related  to  the  objective  function  of  the  model.  The  best 
solutions  are  reserved  for  crossover  and  inferior  solutions  are  eliminated. 


Crossover  of  the  best  solutions  creates  chromosomes  to  replace  the  ones  that  were  eliminated  from  the 
initial  population.  Two  chromosomes  (parents)  are  chosen  for  crossover,  which  results  in  two  new 
chromosomes  (children).  Consider  a  crossover  example  involving  the  chromosomes  shown  in  Figure  4. 
Let  the  3rd  gene  denote  the  position  of  the  crossover  and  suppose  the  length  of  the  crossover  is  3.  As 
shown  in  Figure  5,  the  first  two  genes  from  parent  one  become  the  first  two  genes  in  child  one,  and  the 
last  three  genes  from  parent  two  become  the  last  three  genes  in  child  one.  The  opposite  is  also  true.  That 
is,  the  first  two  genes  from  parent  two  become  the  first  two  genes  in  child  two,  and  the  last  three  genes 
from  parent  one  become  the  last  three  genes  in  child  two.  The  crossover  routine  and  resulting  child 
chromosomes  are  shown  in  Figure  5. 


2 

1 

0 

0 

1 

1 

0 

2 

1 

1 

Figure  4:  Parent  Chromosomes  Chosen  for  Crossover 


Figure  5:  Crossover  Routine  and  Resulting  Child  Chromosomes 
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Occasionally,  a  GA  may  become  “stuck”  on  a  local  optimal  solution.  To  force  the  algorithm  to  leave  local 
optima  in  search  of  better  solutions,  some  of  the  new  chromosomes  undergo  mutation.  Mutation  is 
performed  by  selecting  a  gene  within  the  chromosome  and  assigning  a  new  value  for  that  gene.  New 
chromosomes  created  with  crossover  and  mutation  are  merged  with  the  parent  chromosomes  to  create  a 
new  generation.  The  GA  is  executed  for  a  specified  number  of  generations. 

In  the  GA  used  to  evaluate  model  P,  the  individual  genes  within  a  chromosome  represent  the  number  of 
components  to  repair  in  subsystem  ij  (that  is,  the  decision  variable,  d„).  Therefore,  each  chromosome 
contains  q  *  m  genes.  An  initial  population  of  candidate  solutions  (chromosomes)  is  created.  The  size  of 
the  population  is  an  input  parameter  specified  by  the  model  user.  To  initialize  the  population,  a  specified 
number  of  chromosomes  are  created  such  that  the  value  of  each  gene  is  a  discrete  uniform  random 
variable  over  the  integers  {0,  1,  ...  ,  «/;}.  The  population  is  initialized  exactly  once  for  each  scenario 
evaluated. 

The  initial  population  represents  the  first  generation.  The  feasibility  of  each  chromosome  in  the 
generation  is  checked  against  model  constraints.  If  a  chromosome  is  found  to  be  infeasible,  it  undergoes 
mutation.  That  is,  a  gene  is  chosen  at  random,  and  the  new  value  of  the  gene  is  initialized.  The 
chromosome  continues  to  undergo  mutation  until  it  is  feasible.  Once  a  generation  of  feasible  solutions 
exists,  each  chromosome  is  ranked  according  to  the  fitness  function.  In  this  case,  the  fitness  function  is 
the  overall  reliability  resulting  from  each  chromosome.  The  best  solutions,  or  top-half  of  the  ranked 
solutions,  are  reserved  for  crossover,  and  the  inferior  solutions  are  eliminated. 

In  the  crossover  routine,  two  parent  chromosomes  are  chosen  at  random.  The  position  and  length  of  the 
crossover  are  also  chosen  randomly.  The  crossover  is  performed  as  shown  in  Figure  4  and  Figure  5.  For 
each  chromosome  created  in  the  crossover  routine,  the  probability  of  mutation  is  defined  to  be  0.1.  That 
is,  each  chromosome  has  a  10%  chance  of  undergoing  mutation.  Once  crossover  and  mutation  are 
complete,  the  parents  and  children  are  merged  to  create  a  new  generation. 

The  new  generation  is  ranked  according  to  overall  reliability.  The  best  solutions  are  subjected  to 
crossover,  and  chromosomes  resulting  from  crossover  are  subjected  to  mutation.  The  parents  and  children 
are  then  merged,  forming  a  new  generation.  This  process  is  repeated  for  a  user-specified  number  of 
generations.  This  logic  is  summarized  in  Figure  6.  This  approach  uses  an  application  developed  in  Visual 
Basic®  and  programmed  as  a  macro  within  an  Excel  spreadsheet.  Instructions  on  using  both  the  total 
enumeration  strategy  and  the  GA  for  solving  model  P  are  presented  in  Appendix  A. 
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4  The  Static  Model 


The  scenario  captured  by  the  optimization  model  P  addresses  selective  maintenance  for  a  set  of  systems, 
but  it  does  not  capture  the  assignment  of  systems  to  missions.  In  the  scenario  address  by  model  P, 
assignment  is  not  relevant,  because  all  missions  in  the  upcoming  set  are  identical.  In  this  section,  we  add 
assignment  to  the  model  by  considering  the  scenario  in  which  the  missions  in  the  upcoming  set  are  not 
identical. 

4. 1  Hypothetical  Set  of  Systems 

Consider  a  set  of  q  independent  and  identical  systems.  Each  system  is  comprised  of  m  independent 
subsystems  connected  in  series,  and  each  subsystem,  say  subsystem  j,  contains  nt  independent  and 
identical  copies  of  a  CFR  component  connected  in  parallel.  A  graphical  representation  (reliability  block 
diagram)  of  an  example  set  of  systems  is  shown  in  Figure  7.  Note  that  each  component  in  Figure  7  is 
labeled  j,  k where  i  denotes  the  system  number,  j  denotes  the  subsystem  number  and  k  denotes  the 
component  number. 

The  set  of  systems  is  required  to  perform  sequential  sets  of  missions  that  originate  from  and  return  to  a 
common  base  of  operation.  Each  set  of  missions  contains  q  independent  missions  with  common  start 
times;  however,  the  missions  have  different  durations.  At  any  time,  a  specific  component,  subsystem  or 
system  is  in  one  of  two  states:  functioning  (1)  or  failed  (0).  Components  only  fail  during  missions  and 
failed  components  can  be  repaired  only  during  the  breaks  between  sets  of  missions.  Note  that  all 
maintenance  is  performed  at  the  base. 
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Figure  7:  Example  Set  of  Systems 

4.2  System  Performance 

Suppose  that  a  set  of  missions  has  just  ended  and  all  systems  have  returned  to  their  base  of  operation  and 
maintenance.  Note  that  component  ijk  refers  to  component  k  in  subsystem  j  of  system  /,  and  note  that 
subsystem  ij  refers  to  subsystem  j  in  system  i.  Let  r»  denote  the  status  of  component  ijk,  let  Yy  denote  the 
status  of  subsystem  ij,  and  let  7,  denote  the  status  of  system  i.  Since  each  subsystem  is  a  parallel 
arrangement  of  its  components, 

k= 1  k=\ 

Equation  4.1 

Since  each  system  is  a  series  arrangement  of  its  subsystems, 

m 

7=1 

Equation  4.2 

Note  that  if  Yi}  ~  0  for  some  j  =■  1,2,  . . w,  then  Yt  —  0  and  system  i  did  not  successfully  complete  the  last 
mission. 
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Prior  to  the  next  set  of  missions,  maintenance  (repair)  can  be  performed  on  some  or  all  of  the  failed 
components.  Let  Xijk  denote  the  status  of  component  ijk  at  the  start  of  the  next  mission,  let  Xy  denote  the 
status  of  subsystem  ij  at  the  start  of  the  next  mission,  and  let  X,  denote  the  status  of  system  i  at  the  start  of 
the  next  mission.  Note  that 


Equation  4.3 


and 


m 

>1 

Equation  4.4 

Furthermore,  note  that  if  Yyk  —  0  and  Xljk  =  1 ,  then  component  ijk  was  repaired  during  the  break  between 
missions.  Also,  note  that  if Xy  =  0  for  some  j  =  1,2,  ...,  m,  then  X,  =  0  and  system  i  is  unable  to  undertake 
its  next  mission. 

Clearly,  the  Xijk  values  have  a  direct  impact  on  the  status  of  components,  subsystems  and  systems  at  the 
end  of  the  upcoming  mission.  Let  Yl  denote  the  status  of  component  ijk  at  the  end  of  the  next  mission, 

let  Y?  denote  the  status  of  subsystem  ij  at  the  end  of  the  next  mission,  and  let  Y*  denote  the  status  of 
system  i  at  the  end  of  the  next  mission.  Note  that 


and 


k~\ 

Equation  4.5 


ffl 

r,*=nv 

7=1 

Equation  4.6 
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Also,  note  that  if  X,jk  =  0,  then  Yl  =  0 .  If  X,,k  =  1,  then  Y*k  is  a  random  variable.  Assuming  A'i  =  1,  then 
both  YjJ for  all  j  =  1,  2,  m  and Y+ ,  i  =  1,  2,  q  are  random  variables.  Furthermore,  note  that  if 

Y*  =0  for  some  j  -  1,  2,  ...,  m,  then  Yt+  =0  and  system  i  did  not  successfully  complete  the  next 
mission. 


Since  the  upcoming  missions  are  not  identical,  the  length  of  the  mission  must  be  incorporated  into  the 
component  reliability  calculation.  In  other  words,  we  must  adjust  the  component  reliability  values 
depending  upon  the  mission  assignment.  Let  A ,•  denote  the  failure  rate  of  a  type  j  component,  and  let  prj 
denote  the  probability  that  a  functioning  type  j  component  completes  its  next  mission  if  it  is  assigned 
mission  t.  Let  L  denote  the  length  of  mission  f,  f  =  1,  2,  ...  ,  q.  Then 


Equation  4.7 

This  parameter  quantifies  the  difficulty  of  each  mission  and  the  resulting  differences  between  missions. 
Here,  we  consider  mission  length  to  be  synonymous  with  mission  difficulty,  but  other  measures  of 
difficulty  could  be  applied  by  adjusting  the  values  of  /?,-/. 

Next,  we  define  conditional  component  reliability,  ry  as  the  probability  that  a  functioning  component  in 
subsystem  ij  completes  its  next  mission 


rf-Z**r 

/■=i 

Equation  4.8 

where  g,y  denotes  the  mission  assignment  variable  and  is  given  by 


fl  if  system  i  is  assigned  to  mission  f ' 
[0  otherwise 

Equation  4.9 
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Note  that  conditional  component  reliability  is  only  indexed  on  i  and  j  because  the  components  within  a 
subsystem  are  identical.  The  resulting  measures  of  component  and  subsystem  reliability  are  given  by 


^ijk  ^"ij^ ijk 

Equation  4.10 

Equation  4.11 


where 


b  =  n 

a  i 


■aii+du 


Equation  4.12 

Finally,  the  reliability  of  system  i  is  defined  to  be  the  probability  that  system  /  is  functioning  at  the  end  of 
the  next  mission  and  given  by 


m 

Equation  4.13 

4.3  Formulation  of  the  Static  Model 

Recall  that  we  consider  the  case  in  which  a  set  of  missions  has  just  ended  and  all  systems  have  returned  to 
their  base  of  operation  and  maintenance.  The  base  possesses  the  technological  capability  to  repair  any 
failed  component  and  ideally,  all  failed  components  would  be  repaired  prior  to  the  beginning  of  the  next 
mission;  however,  maintenance  resource  limitations  may  prevent  the  repair  of  all  failed  components. 

Let  aij  denote  the  number  of  failed  components  in  subsystem  ij\  and  note  that  a is  given  by 

*-l 

Equation  4.14 
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Furthermore,  the  total  number  of  failed  components  of  type  j  is  given  by 


a-j=Tau 

/= 1 

Equation  4.15 

Each  repair  consumes  a  fixed  amount  of  each  of  s  maintenance  resources  (labor,  equipment,  and  spare 
parts,  for  example.).  Let  aji  denote  the  amount  of  resource  /  consumed  by  repairing  a  component  of  type y, 
and  suppose  /?/  denotes  the  amount  of  resource  l  available  during  a  single  break.  If 

m 

.H 

Equation  4.16 

for  all  /  =  1,  2,  ...  ,  s,  then  all  failed  components  may  be  repaired  prior  to  the  next  mission.  Otherwise,  a 
method  is  needed  to  decide  which  failed  components  should  be  repaired  and  which  components  must 
remain  in  a  failed  condition.  This  is  the  selective  maintenance  scenario  that  we  consider. 

Let  djj  denote  the  number  of  failed  components  in  subsystem  ij  to  repair  prior  to  the  beginning  of  the  next 
mission.  Then 


*=t 


Equation  4.17 

The  total  number  of  failed  components  of  type  j  repaired  prior  to  the  next  mission  is  given  by 

/=1 

Equation  4.18 

To  incorporate  assignment  into  the  model,  the  mission  assignment  variable,  gif%  is  defined.  To  address  the 
selective  maintenance  issue,  we  formulate  a  nonlinear,  discrete  optimization  model  where  the  dy  and  g & 
values  are  the  decision  variables. 
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The  first  constraint  on  the  dy  decision  variables  is  that  they  be  integer-valued.  Second,  the  number  of 
repairs  is  limited  to  the  number  of  failed  components,  that  is, 

O^dy^Cly 

Equation  4.19 

Third,  the  number  of  repairs  is  limited  by  the  available  maintenance  resources,  specifically, 

rn 

Equation  4.20 


for  all  /=  1, 2,  5. 

The  first  constraint  on  the  gfr  decision  variables  is  that  they  are  binary  variables.  Second,  a  system  must 
be  assigned  exactly  one  mission,  that  is, 


]L£*'=1 

Equation  4.21 


Third,  each  mission  must  be  assigned  to  exactly  one  system,  specifically, 

</ 

YuSrr 

i"=i 

Equation  4.22 
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The  objective  in  choosing  the  values  for  the  decision  variables  is  to  maximize  the  probability  that  all 
upcoming  missions  are  completed  successfully.  This  probability  is  referred  to  as  overall  reliability  and  is 
given  by 

^iWnrM-'v)** 

/=]  /=!  ./-I 


Equation  4.23 


The  full  formulation  of  the  optimization  model  for  this  scenario  is  as  follows: 


SP:  Maximize  ^  =  Fini”  ^  ^  Sir 


<L _ m  (  ±_ 

> 

<=1  J= 1  V  <"= 


subject  to 


q  m 

-  Pi 


M  ./=! 


/=  1, 2,  ,  5 


±g,=  1 


I>,=1  ' 

/"=! 

gw ,  binary  i=\,2,...,q 


0  <  d(j  <  a(j ,  integer  i=\,2,...,q 

j=  1,2, ... ,  m 

Equation  4.24 
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4.4  Solution  Procedure 


Due  to  the  complexity  of  incorporating  both  selective  maintenance  and  mission  assignment  into  the 
model,  we  prefer  to  use  a  search-based  heuristic  such  as  GAs  rather  than  total  enumeration.  GAs  are 
computationally  intense  search  methods  originally  applied  in  the  area  of  artificial  intelligence.  Over  the 
past  decade,  increases  in  computer  speeds  have  made  this  directed  trial-and-error  method  increasingly 
popular.  GAs  operate  on  the  principle  of  “survival  of  the  fittest.”  In  general,  GA  keeps  track  of  a  fixed 
population  of  candidate  solutions  called  chromosomes.  Each  element,  or  gene,  within  a  chromosome 
represents  some  portion  of  the  solution.  The  solutions  are  then  ranked  using  a  fitness  function  related  to 
the  objective  function  of  the  model.  The  best  solutions  are  reserved  for  crossover,  and  inferior  solutions 
are  eliminated. 

Crossover  of  the  best  solutions  creates  chromosomes  to  replace  the  ones  that  were  eliminated  from  the 
initial  population.  Two  chromosomes  (parents)  are  chosen  for  crossover,  which  results  in  two  new 
chromosomes  (children).  This  combination  of  the  parents  and  children  is  referred  to  as  a  generation. 

Occasionally,  a  GA  may  become  “stuck”  on  a  local  optimal  solution.  To  force  the  algorithm  to  leave  local 
optima  in  search  of  better  solutions,  some  of  the  new  chromosomes  undergo  mutation.  Mutation  is 
performed  by  selecting  a  gene  within  the  chromosome  and  assigning  a  new  value  for  that  gene.  New 
chromosomes  created  with  crossover  and  mutation  are  merged  with  the  parent  chromosomes  to  create  a 
new  generation.  The  GA  repeats  itself  for  a  specified  number  of  generations. 

In  the  GA  used  to  evaluate  model  SP,  the  individual  genes  within  a  chromosome  represent  not  only  the 
mission  assignment  (that  is,  the  decision  variable  g„)  but  also  the  number  of  components  to  repair  in 
subsystem  ij  (specifically,  the  decision  variable  du).  Therefore,  each  chromosome  contains  q  +  q*  m 
genes.  An  initial  population  of  candidate  solutions  (chromosomes)  is  created.  The  size  of  the  population  is 
an  input  parameter  specified  by  the  model  user.  To  initialize  the  population,  a  specified  number  of 
chromosomes  are  created  such  that  the  value  of  the  first  q  genes  are  discrete  uniform  random  variables 
over  the  integers  {1,2,  ...  ,  q).  The  values  of  these  genes  are  generated  by  randomly  selecting  one  of  the 
q\  feasible  mission  assignments.  The  value  of  the  remaining  q*m  genes  represent  the  selective 
maintenance  decisions.  The  value  of  these  genes  are  generated  as  discrete  uniform  random  variables  over 
the  integers  {0,  1,  ...  ,  a,/} .  The  population  is  initialized  exactly  once  for  each  scenario  evaluated. 
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The  initial  population  represents  the  first  generation.  The  feasibility  of  each  chromosome  in  the 
generation  is  checked  against  model  constraints.  If  a  chromosome  is  found  to  be  infeasible,  two  checks 
are  made.  First,  the  model  checks  to  see  if  any  duplicate  mission  assignments  occur  in  the  first  q  genes.  If 
a  duplicate  mission  assignment  is  found,  one  of  the  duplicate  genes  is  chosen  at  random  and  initialized 
again.  The  chromosome  is  checked  again  for  duplicate  missions,  and  the  above  process  is  repeated  until 
no  duplicate  mission  assignments  exist.  Second,  the  values  of  the  remaining  q*m  chromosomes  may 
violate  model  constrains.  If  this  is  the  case,  the  chromosome  continues  to  undergo  mutation  until  it  is 
feasible.  That  is,  one  gene  is  chosen  at  random  and  initialized  again.  Once  a  generation  of  feasible 
solutions  exists,  each  chromosome  is  ranked  according  to  a  fitness  function.  In  this  case,  the  fitness 
function  is  the  overall  reliability  resulting  from  each  chromosome.  The  best  solutions  are  reserved  for 
crossover,  and  the  inferior  solutions  are  eliminated. 

In  the  crossover  routine,  two  parent  chromosomes  are  chosen  at  random.  In  this  scenario,  crossover 
occurs  in  two  steps.  First,  crossover  is  performed  on  the  first  q  genes  representing  the  mission  assignment. 
After  this  portion  of  crossover  is  complete,  the  genes  are  checked  for  duplicate  assignments.  Considering 
the  chromosomes  shown  in  Figure  §,  suppose  the  position  of  the  crossover  is  1,  and  the  length  of  the 
crossover  is  2.  The  chromosomes  after  the  first  crossover  step  are  shown  in  Figure  9 . 


Figure  8:  Parent  Chromosomes  Chosen  for  Crossover 
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Figure  9:  Step  One  of  Crossover  Routine 


Notice  that  in  each  of  the  new  chromosomes,  a  mission  assignment  is  duplicated.  Specifically,  the 
chromosome  on  the  right  has  two  systems  assigned  to  mission  three  and  none  to  mission  one,  and  the 
chromosome  on  the  left  has  two  systems  assigned  to  mission  one  and  none  to  mission  three.  To  correct 
this  situation,  a  duplicate  assignment  gene  is  chosen  from  each  chromosome.  The  value  of  that  gene  is 
initialized  again.  The  chromosome  is  checked  again  for  duplicate  missions,  and  the  above  process  is 
repeated  until  no  duplicate  mission  assignments  exist. 

Next,  the  remaining  q*m  genes  undergo  crossover.  The  position  and  length  of  the  crossover  are  also 
chosen  randomly.  Again,  considering  the  chromosomes  shown  in  Figure  8 ,  suppose  the  position  of  the 
crossover  is  gene  2  (in  the  second  section  of  the  chromosome)  and  the  length  of  the  crossover  is  3.  Figure 
10  shows  the  second  step  in  the  crossover  routine. 
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Figure  10:  Step  Two  of  Crossover  Routine 


For  each  chromosome  created  in  the  crossover  routine,  the  probability  of  mutation  is  defined  to  be  0.1. 
That  is,  each  chromosome  has  a  10%  chance  of  undergoing  mutation.  Once  crossover  and  mutation  are 
complete,  the  parents  and  children  are  merged  to  create  a  new  generation. 

The  new  generation  is  ranked  according  to  overall  reliability.  The  best  solutions  are  subjected  to 
crossover,  and  chromosomes  resulting  from  crossover  are  subjected  to  mutation.  The  parents  and  children 
are  then  merged,  forming  a  new  generation.  This  process  is  repeated  for  a  user-specified  number  of 
generations. 

As  a  numerical  example,  consider  the  system  presented  in  Figure  77,  having  q  -  4  and  m  =  3.  Two  ( s  ~  2) 
limited  maintenance  resources  limit  repair  activities.  Note  that  J3\  =  16  and  /%  =  14.  The  next  mission 
durations  (in  hours)  are  t\  =  4,  t2  =  1 .75,  t3  =  2,  t4  -  0.5.  The  remaining  parameters  for  the  set  of  systems 
are  given  in  Table  4.1.  To  repair  all  failed  components  requires  30  units  of  resource  1  and  21  units  of 
resource  2.  Therefore,  we  can  use  our  selective  maintenance  model  to  determine  the  mission  assignment 
for  each  system  and  which  components  to  repair  to  maximize  the  overall  reliability.  We  solve  this 
scenario  using  the  GA.  This  approach  uses  an  application  developed  in  Visual  Basic®  and  programmed  as 
a  macro  within  a  Microsoft®  Excel  spreadsheet.  Instructions  for  using  the  spreadsheet  appear  in  Appendix 
B.  Figure  12  shows  the  system  status  after  maintenance.  The  overall  reliability  of  the  set  of  systems  with 
the  proposed  repairs  and  mission  assignment  is  0.5347.  Given  this  mission  assignment,  overall  reliability 
would  be  0.6039  if  sufficient  maintenance  resources  were  available  to  repair  all  failed  components. 
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Table  4.1:  Model  Parameters  for  Example  Set  of  Systems 


system  1 


mission  4 


system  2 


mission  1  system  3 


component  not  repaired  during  break 


component  repaired  during  break 


Figure  12:  Mission  Assignment  and  Status  after  Maintenance 
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5  The  Dynamic  Model 

The  scenario  captured  by  the  optimization  model  P  addresses  selective  maintenance  for  a  set  of  systems 
but  not  assignment.  Assignment  is  not  addressed  because  all  missions  in  the  upcoming  set  are  identical.  In 
Section  4,  we  add  assignment  to  the  model  by  considering  the  scenario  in  which  the  missions  in  the 
upcoming  set  are  not  identical;  however,  all  missions  within  a  set  have  common  start  times.  In  this 
scenario,  missions  start  and  end  at  different  times,  and  maintenance  and  scheduling  decisions  are  made 
over  a  series  of  time  “buckets.” 

5.1  Hypothetical  Set  of  Systems 

Consider  a  set  of  q  independent  and  identical  systems.  Each  system  is  comprised  of  m  independent 
subsystems  connected  in  series,  and  each  subsystem,  say  subsystem  j,  contains  nj  independent  and 
identical  copies  of  a  CFR  component  connected  in  parallel.  A  graphical  representation  (reliability  block 
diagram)  of  an  example  set  of  systems  is  shown  in  Figure  13.  Note  that  each  component  in  Figure  13  is 
labeled  j,  k ,”  where  i  denotes  the  system  number,  /  denotes  the  subsystem  number  and  k  denotes  the 
component  number. 


Figure  13:  Example  Set  of  Systems 

The  set  of  systems  is  required  to  perform  sequential  sets  of  missions  that  originate  from  and  return  to  a 
common  base  of  operation.  Each  set  of  missions  contains  q  independent  missions.  Missions  have  varying 
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start  times  and  durations.  Additionally,  systems  do  not  necessarily  complete  their  previous  missions  at  the 
same  time.  At  any  time,  a  specific  component,  subsystem  or  system  is  in  one  of  two  states:  functioning 
(1)  or  failed  (0).  Components  only  fail  during  missions,  and  failed  components  can  be  repaired  only 
during  the  breaks  between  sets  of  missions.  Note  that  all  maintenance  is  performed  at  the  base. 

5.2  System  Performance 

Suppose  that  system  i  returns  to  its  base  of  operation  and  maintenance  from  its  current  mission  at  the 
beginning  of  time  bucket  RT{.  We  define  “now”  as  the  beginning  of  time  bucket  1,  and  at  least  one  system 
has  RTj  <  1 .  During  the  planning  horizon,  all  systems  not  initially  at  the  base  will  return  from  their  current 
mission. 

Note  that  component  ijk  refers  to  component  k  in  subsystem  j  of  system  /,  and  note  that  subsystem  ij 
refers  to  subsystem  j  in  system  j.  Let  denote  the  status  of  component  ijk,  let  Yy  denote  the  status  of 
subsystem  ij,  and  let  7,  denote  the  status  of  system  /.  Since  each  subsystem  is  a  parallel  arrangement  of  its 
components, 

k=\  k=] 


Equation  5.1 

Since  each  system  is  a  series  arrangement  of  its  subsystems, 

m 

rl=nr» 

j*  i 

Equation  5.2 

Note  that  if  Yy  =  0  for  some  j  =  1 , 2, . . .,  m,  then  7,  =  0  and  system  i  did  not  successfully  complete  the  last 
mission. 

Prior  to  the  system’s  next  assigned  mission,  maintenance  (repair)  can  be  performed  on  some  or  all  of  the 
failed  components.  Let  X,jk  denote  the  status  of  component  ijk  at  the  start  of  the  next  mission,  let  Xy 
denote  the  status  of  subsystem  ij  at  the  start  of  the  next  mission,  and  let  X,  denote  the  status  of  system  i  at 
the  start  of  the  next  mission.  Note  that 
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x,=Ux 


<jk 


k-\ 


Equation  5.3 

and 

m 

*,=n*. 

>1 


Equation  5.4 

Furthermore,  note  that  if  V=  0  and  Xyk  =  1 ,  then  component  ijk  was  repaired  during  the  break  between 
missions.  Also,  note  that  if  Xy  =  0  for  some  j  =  1,2,  . . .,  m,  then  Xj  =  0  and  system  /  is  unable  to  undertake 
its  next  mission. 

Clearly,  the  A'p  values  have  a  direct  effect  on  the  status  of  components,  subsystems  and  systems  at  the 
end  of  the  upcoming  mission.  Let  Yl  denote  the  status  of  component  ijk  at  the  end  of  the  next  mission, 

let  YjJ  denote  the  status  of  subsystem  ij  at  the  end  of  the  next  mission,  and  let  Yl  denote  the  status  of 
system  i  at  the  end  of  the  next  mission.  Note  that 


Y+ 

!l 


k=\ 


Equation  5.5 

and 

m 

i?=njy 

7=1 

Equation  5.6 


Also,  note  that  if Xyu  =  0,  then  Yl  =  0  .  If XiJk  =  1,  then  Yl  is  a  random  variable.  Assuming^  =  1,  then 
both  l^+for  all  j  =  1,  2,  ...,  m  and}^+,  /  =  1,  2,  ...,  q  are  random  variables.  Furthermore,  note  that  if 
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YjJ  =  0  for  some  j  =  1,  2,  m,  then  Y+  =  0  and  system  i  did  not  successfully  complete  the  next 
mission. 

Each  system  within  the  set  must  be  assigned  exactly  one  upcoming  mission.  Mission  f  starts  at  the 
beginning  of  time  bucket  STp,  and  tf  denotes  the  length  of  mission  i .  Note  that  no  future  mission  ends 
before  the  end  of  the  planning  horizon.  Since  the  upcoming  missions  are  not  identical,  the  length  of  the 
mission  must  be  incorporated  into  the  component  reliability  calculation.  In  other  words,  we  must  adjust 
the  component  reliability  values  depending  upon  the  mission  assignment.  Let  Xj  denote  the  failure  rate  of 
a  type  j  component,  j  =1,2 , ,m,  and  let  p/y  denote  the  probability  that  a  functioning  type  j  component 
completes  its  next  mission  if  it  is  assigned  mission  f,!=  1,2,  ...  ,  q,j  =  1,  2,  ...  ,  m.  Let  tr  denote  the 
length  of  mission  f,  t  =  1,  2, ...  ,  q. 


Equation  5.7 

This  parameter  quantifies  the  difficulty  of  each  mission  and  the  resulting  differences  between  missions. 
Here,  we  consider  mission  length  to  be  synonymous  with  mission  difficulty,  but  other  measures  of 
difficulty  could  be  applied  by  adjusting  the  values  of p,y,  i=  1,2,...,  qj  =1,2 , ,m. 

Next,  we  define  conditional  component  reliability,  ry,  as  the  probability  that  a  functioning  component  in 
subsystem  ij completes  its  next  mission,  /=  1,2,  ...  ,  q,j=  1,2,  ...  ,  m 

'/ 

r,j  PrjSw 

/•= i 

Equation  5.8 

where  glf  denotes  the  mission  assignment  variable  and  is  given  by 

fl  if  system  /  is  assigned  to  mission  /' 
g"  [0  otherwise 

Equation  5.9 
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/  =  1,  2,  ...  ,  q,  f  =  1,  2,  ...  ,  q.  Note  that  conditional  component  reliability  is  only  indexed  on  /  and  j 
because  the  components  within  a  subsystem  are  identical.  The  resulting  measures  of  component  and 
subsystem  reliability  are  given  by 


^ijk  yk 

Equation  5.10 


and 


Equation  5.11 


where 


Equation  5.12 

Finally,  the  reliability  of  system  /  is  defined  to  be  the  probability  that  system  i  is  functioning  at  the  end  of 
the  next  mission  and  given  by 

m 

*,=n  r, 

i=i 

Equation  5.13 

5.3  Formulation  of  the  Dynamic  Model 

Recall  that  we  consider  the  case  in  which  system  i  has  returned  from  its  current  mission  at  the  beginning 
of  time  bucket  RT,.  Its  base  of  operation  and  maintenance  possesses  the  technological  capability  to  repair 
any  failed  component  and  ideally,  all  failed  components  would  be  repaired  prior  to  the  beginning  of  the 
next  mission;  however,  maintenance  resource  limitations  may  prevent  the  repair  of  all  failed  components. 

Let  ai}  denote  the  number  of  failed  components  in  subsystem  ij\  and  note  that  is  given  by 


36 


k=\ 

Equation  5.14 

Furthermore,  the  total  number  of  failed  components  of  type  j  is  given  by 

'/ 

a'j=Hav 

i=i 

Equation  5.15 

Note  that  for  systems  having  RT,  >  1,  the  values  for  the  number  of  failed  components  in  subsystem  ij,  ay 
are  estimates.  Each  repair  consumes  a  fixed  amount  of  each  of  s  maintenance  resources  (labor,  equipment 
and  spare  parts,  for  example).  Let  aj/  denote  the  amount  of  resource  /  consumed  by  repairing  a  component 
of  type  j,  and  suppose  Pi  denotes  the  amount  of  resource  l  available  during  in  each  time  bucket  h.  If 
sufficient  resources  exist,  then  all  failed  components  may  be  repaired  prior  to  the  next  mission. 
Otherwise,  a  method  is  needed  to  decide  which  failed  components  should  be  repaired  and  which 
components  must  remain  in  a  failed  condition.  This  is  the  selective  maintenance  scenario  that  we 
consider. 

Let  diji,  denote  the  number  of  failed  components  in  subsystem  ij  to  repair  in  time  bucket  h.  The  total 
number  (across  all  systems)  of  type  j  components  to  repair  during  time  bucket  h  is  given  by 

</ 

d.jh  =  J]  dijh 

/'= i 

Equation  5.16 

The  total  number  of  components  to  repair  in  subsystem  ij  during  the  planning  horizon  is  given  by 

dy,  —  Yd* 

h= 1 

Equation  5.17 


where 
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-max _  max 

1 


Equation  5.18 

To  incorporate  assignment  into  the  model,  an  assignment  decision  variable,  g,r,  is  defined.  To  address  the 
selective  maintenance  issue,  we  formulate  a  nonlinear  discrete  optimization  model  where  the  dy  and  gir 
values  are  the  decision  variables. 

The  first  constraint  on  the  dy  decision  variables  is  that  they  be  integer-valued.  Second,  the  number  of 
repairs  in  each  time  bucket  must  be  greater  than  or  equal  to  zero.  Third,  the  number  of  repairs  is  limited 
by  the  available  maintenance  resources,  such  that 


'Lajid-j  *  Pi 


M 

Equation  5.19 


for  all  /=  1,2, ...,  s. 

The  first  constraint  on  the  gy  decision  variables  is  that  they  be  binary.  Second,  a  system  must  be  assigned 
exactly  one  mission,  specifically, 


i>,= i 

M 

Equation  5.20 

Third,  each  mission  must  be  assigned  to  exactly  one  system,  that  is. 


Equation  5.21 
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Additional  model  constraints  are  necessary  to  restrict  the  time  buckets  in  which  maintenance  actions  and 
mission  assignments  may  be  performed  on  a  particular  system.  First,  maintenance  cannot  be  performed  if 
a  system  has  not  returned  from  its  current  mission.  That  is, 

q  m  RTj-\ 

II  I4»  =° 

/=i  7=1  h= l 

Equation  5.22 

for  all  h  =  1,2,  ...,  RT/  -  1 .  Second,  a  system  cannot  be  assigned  a  mission  if  it  has  not  returned  from  its 
previous  mission.  This  is  given  by 

g„{RT,  -  ST, -1)<0 


Equation  5.23 

for  all  /  =  1,2,  . . .,  q  and  i  =  1, 2, ...,  q.  Finally,  maintenance  cannot  be  performed  if  a  system  has  started 
its  next  mission,  specifically, 

*,(A-Srf)iw£0 

Equation  5.24 


for  all  i—  1, 2, ...,  q,  f  =  1,  2,  ...,  q ,  and  h  =  1,2,...,  ST” 


The  objective  in  choosing  values  for  the  decision  variables  is  to  maximize  the  probability  that  all 
upcoming  missions  are  completed  successfully.  This  probability  is  referred  to  as  overall  reliability  and  is 
given  by 

/=l  ,=l  7=1 

Equation  5.25 

The  full  formulation  of  the  optimization  model  for  this  scenario  is  as  follows: 
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gw*  binary 


*=1.2, 


1.2, 


5.4  Solution  Procedure 

Due  to  the  complexity  of  incorporating  both  selective  maintenance  and  a  more  dynamic  mission  profile 
into  the  model,  we  prefer  to  use  a  search-based  heuristic  such  as  GAs  rather  than  total  enumeration.  GAs 
are  computationally  intense  search  methods  originally  applied  in  the  area  of  artificial  intelligence.  Over 
the  past  decade,  increases  in  computer  speeds  have  made  this  directed  trial-and-error  method  increasingly 
popular.  GAs  operate  on  the  principle  of  “survival  of  the  fittest.”  In  general,  GA  keeps  track  of  a  fixed 
population  of  candidate  solutions  called  chromosomes.  Each  element,  or  gene,  within  a  chromosome 
represents  some  portion  of  the  solution.  The  solutions  are  then  ranked  using  a  fitness  function  related  to 
the  objective  function  of  the  model.  The  best  solutions  are  reserved  for  crossover  and  inferior  solutions 
are  eliminated. 

Crossover  of  the  best  solutions  creates  chromosomes  to  replace  the  ones  that  were  eliminated  from  the 
initial  population.  Two  chromosomes  (parents)  are  chosen  for  crossover,  which  results  in  two  new 
chromosomes  (children).  This  combination  of  the  parents  and  children  is  referred  to  as  a  generation. 

Occasionally,  a  GA  may  become  “stuck”  on  a  local  optimal  solution.  To  force  the  algorithm  to  leave  local 
optima  in  search  of  better  solutions,  some  of  the  new  chromosomes  undergo  mutation.  Mutation  is 
performed  by  selecting  a  gene  within  the  chromosome  and  assigning  a  new  value  for  that  gene.  New 
chromosomes  created  with  crossover  and  mutation  are  merged  with  the  parent  chromosomes  to  create  a 
new  generation.  The  GA  repeats  itself  for  a  specified  number  of  generations. 

In  the  GA  used  to  evaluate  model  DSP,  the  individual  genes  within  a  chromosome  represent  not  only  the 
mission  assignment  (explicitly,  the  decision  variable  gw),  but  also  the  number  of  components  to  repair  in 
subsystem  ij  during  each  time  “bucket”  (i.e.,  the  decision  variable  d#,),  Therefore,  each  chromosome 
contains  q  +  q*  m*h  genes.  An  initial  population  of  candidate  solutions  (chromosomes)  is  created.  The 
size  of  the  population  is  an  input  parameter  specified  by  the  model  user.  To  initialize  the  population,  a 
specified  number  of  chromosomes  are  created  such  that  the  value  of  the  first  q  genes  are  discrete  uniform 
random  variables  over  the  integers  {1,2,  ...  ,  q).  The  values  of  these  genes  are  generated  by  randomly 
selecting  one  of  the  q\  feasible  mission  assignments.  The  value  of  the  remaining  q*m*h  genes  represent 
the  selective  maintenance  decisions.  The  value  of  these  genes  are  generated  as  discrete  uniform  random 
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variables  over  the  integers  {0,  1,  ...  ,  ay}.  The  population  is  initialized  exactly  once  for  each  scenario 
evaluated. 

The  initial  population  represents  the  first  generation.  The  feasibility  of  each  chromosome  in  the 
generation  is  checked  against  model  constraints.  If  a  chromosome  is  found  to  be  infeasible,  two  checks 
are  made.  First,  the  model  checks  to  see  if  any  duplicate  mission  assignments  occur  in  the  first  q  genes.  If 
a  duplicate  mission  assignment  is  found,  one  of  the  duplicate  genes  is  chosen  at  random  and  initialized 
again.  The  chromosome  is  checked  again  for  duplicate  missions,  and  the  above  process  is  repeated  until 
no  duplicate  mission  assignments  exist.  Second,  the  values  of  the  remaining  q*m*h  chromosomes  may 
violate  model  constraints.  If  this  is  the  case,  the  chromosome  continues  to  undergo  mutation  until  it  is 
feasible.  In  this  case,  we  define  procedures  for  mutation  that  attempt  to  find  a  feasible  solution.  That  is, 
one  gene  whose  value  is  greater  than  or  equal  to  1  is  chosen  at  random  and  the  value  is  decremented  by 
one.  Once  a  generation  of  feasible  solutions  exists,  each  chromosome  is  ranked  according  to  a  fitness 
function.  In  this  case,  the  fitness  function  is  the  overall  reliability  resulting  from  each  chromosome.  The 
best  solutions  are  reserved  for  crossover,  and  the  inferior  solutions  are  eliminated. 

In  the  crossover  routine,  two  parent  chromosomes  are  chosen  at  random.  In  this  scenario,  crossover 
occurs  in  two  steps.  First,  crossover  is  performed  on  the  first  q  genes  representing  the  mission  assignment. 
After  this  portion  of  crossover  is  complete,  the  genes  are  checked  for  duplicate  assignments.  Considering 
the  chromosomes  shown  in  Figure  14,  suppose  the  position  of  the  crossover  is  1 ,  and  the  length  of  the 
crossover  is  2.  The  chromosomes  after  the  first  crossover  step  are  shown  in  Figure  15. 


Figure  14:  Parent  Chromosomes  Chosen  for  Crossover 
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Figure  15:  Step  One  of  Crossover  Routine 


Notice  that  in  each  of  the  new  chromosomes,  a  mission  assignment  is  duplicated.  Specifically,  the 
chromosome  on  the  right  has  two  systems  assigned  to  mission  three  and  none  to'  mission  one,  and  the 
chromosome  on  the  left  has  two  systems  assigned  to  mission  one  and  none  to  mission  three.  To  correct 
this  situation,  a  duplicate  assignment  gene  is  chosen  from  each  chromosome.  The  value  of  that  gene  is 
initialized  again.  The  chromosome  is  checked  again  for  duplicate  missions,  and  the  above  process  is 
repeated  until  no  duplicate  mission  assignments  exist. 

Next,  the  remaining  q*m*h  genes  undergo  crossover.  The  position  and  length  of  the  crossover  are  also 
chosen  randomly.  Again,  consider  the  chromosomes  shown  in  Figure  14.  Suppose  the  position  of  the 
crossover  is  gene  2  (in  the  second  section  of  the  chromosome)  and  the  length  of  the  crossover  is  3.  Figure 
16  shows  the  second  step  in  the  crossover  routine. 
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Figure  16:  Step  Two  of  Crossover  Routine 


For  each  chromosome  created  in  the  crossover  routine,  the  probability  of  mutation  is  defined  to  be  0.1. 
That  is,  each  chromosome  has  a  10%  chance  of  undergoing  mutation.  Once  crossover  and  mutation  are 
complete,  the  parents  and  children  are  merged  to  create  a  new  generation. 

The  new  generation  is  ranked  according  to  overall  reliability.  The  best  solutions  are  subjected  to 
crossover,  and  chromosomes  resulting  from  crossover  are  subjected  to  mutation.  The  parents  and  children 
are  then  merged  forming  a  new  generation.  This  process  is  repeated  for  a  user-specified  number  of 
generations. 

As  a  numerical  example,  consider  the  dynamic  mission  profile  shown  in  Figure  17.  In  this  scenario,  we 
consider  the  set  of  systems  shown  in  Figure  13  have  q  —  4  and  m  -  3.  We  consider  s  =  2  limited 
maintenance  resources,  specifically,  =  8  and  =  6.  The  next  mission  durations  (in  hours)  are  t\  =  4,  t2 
=  1.75,  —  2,  /,}  =  0.5,  and  the  remaining  system  parameters  are  shown  in  Table  5.1.  Note  that  system  1 

returns  from  its  previous  mission  at  the  beginning  of  time  bucket  1 .  Systems  2  and  4  will  return  at  the 
beginning  of  time  bucket  3,  and  system  3  will  return  at  the  beginning  of  time  bucket  4.  Further,  note  that 
upcoming  missions  1  and  2  depart  at  the  beginning  of  time  bucket  4,  mission  3  at  the  beginning  of  time 
bucket  5,  and  mission  4  at  the  beginning  of  time  bucket  6.  This  implies  that  STnax  =  6. 
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opportunities  for 


time  buckets 


Figure  17:  Example  Dynamic  Mission  Profile 


Table  5.1:  Model  Parameters  for  Example  Set  of  Systems 


Subsystem 

/Lj{in  failures  per  Hour) 

■  •• 

aj2 

1 

0.10 

2 

3 

2 

0.16 

3 

1 

3 

0.06 

2 

i 

2 

At  the  beginning  of  time  bucket  1,  “now,”  system  1  returns  from  its  mission.  This  system  has  three  failed 
components  as  showm  in  Figure  18.  At  this  point,  the  number  of  failed  components  in  the  remaining 
systems  are  just  estimates. 
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failed  components  by  subsystem 


Figure  18:  Dynamic  Mission  Profile  with  Component  Failures 


We  can  use  our  selective  maintenance  model  to  determine  the  mission  assignment  for  each  system  and 
which  components  to  repair  in  each  time  bucket  to  maximize  the  overall  reliability.  We  solve  this 
scenario  using  the  GA  described  in  Section  5.4.  This  approach  uses  an  application  developed  in  Visual 
Basic®  and  programmed  as  a  macro  within  an  Excel  spreadsheet.  Instructions  for  using  the  spreadsheet 
appear  in  Appendix  C.  The  model  outputs  are  shown  in  Figure  19;  specifically,  system  1  is  assigned  to 
upcoming  mission  2.  It  will  have  a  component  in  subsystem  1  repaired  in  time  bucket  1  and  a  component 
in  subsystem  3  repaired  in  time  bucket  2.  System  2  is  assigned  to  upcoming  mission  3.  It  will  have  a 
component  in  subsystem  1  repaired  in  time  bucket  3  and  a  component  in  subsystem  2  repaired  in  time 
bucket  4.  System  3  is  assigned  to  upcoming  mission  4.  It  will  have  one  component  in  subsystem  3 
repaired  in  time  bucket  4  and  another  one  in  time  bucket  5.  Finally,  system  4  is  assigned  upcoming 
mission  1  and  will  have  a  component  in  subsystem  2  repaired  in  time  bucket  2.  The  overall  reliability  of 
this  solution  is  0.5225. 
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next  mission  assignments 


in  time  bucket  1 ,  repair  a  type  1 
component  in  system  1 


in  time  bucket  4,  repair  a  type  2 
component  in  system  2  and  a 
type  3  component  in  system  3 


Figure  19:  Solution  for  Numerical  Example 


It  is  important  to  note  that  the  solution  presented  was  run  at  the  beginning  of  time  bucket  1,  and  the 
number  of  failed  components  in  systems  2,  3,  and  4  were  estimates.  At  the  beginning  of  time  bucket  2,  the 
number  of  failed  components  in  systems  2  and  4  will  be  known.  Therefore,  it  may  be  beneficial  to  enter 
the  updated  information  into  the  model  and  run  it  again. 
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6  Future  Opportunities 

Vast  opportunities  for  extending  this  research  exist.  For  example,  all  three  of  the  models  presented  could 
be  extended  to  include  a  variable  cost  of  maintenance.  The  model’s  objective  function  could  be  modified 
to  minimize  the  cost  associated  with  maintenance  actions;  however,  such  a  modification  would  require  a 
penalty  value  associated  with  failed  missions.  Another  extension  could  include  the  ability  to  cancel  a 
mission  if  its  probability  of  success  is  below  some  threshold. 

The  selective  maintenance  models  presented  in  this  research  treat  decision-making  relative  to  a  single 
future  mission.  These  models  could  be  extended  to  capture  the  reliability  and  performance  of  the  set  of 
systems  over  multiple  missions.  If  a  system  is  required  to  perform  a  sequence  of  missions,  then  the 
selective  maintenance  decisions  directly  affect  system  reliability  for  the  next  mission  and  indirectly  affect 
the  system  reliability  for  later  missions.  This  multiple-mission  selective  maintenance  issue  will  be 
addressed  in  Applied  Logistics  Research  Project  BSIT0204.  The  models  presented  in  this  research  will 
serve  as  a  foundation  for  the  upcoming  project. 

Also,  the  models  presented  in  this  research  consider  components,  subsystems,  and  systems  that  are  in  one 
of  two  states:  functioning  or  failed.  It  may  be  more  realistic  to  classify  component  status  using  more  than 
two  discrete  levels.  Additionally,  the  performance  of  a  system  can  be  evaluated  using  several  measures. 
Therefore,  another  potential  extension  to  this  research  includes  considering  multiple  states  of  component 
function  and  multiple  performance  measures  for  a  set  of  systems. 

Currently,  issues  surrounding  cannibalization  are  being  investigated  in  TLI-Military  Logistics  research 
project  MM0206.  Perhaps  in  future  efforts,  selective  maintenance  modeling  could  be  incorporated  into 
cannibalization  decisions. 

Finally,  we  would  like  to  develop  a  USAF-specific  example  to  apply  the  models  developed  in  this 
research  and  the  research  conducted  in  the  other  projects  mentioned  above.  Through  continuous 
collaboration  with  our  partners  at  the  Air  Force  Research  Laboratory  (AFRL),  we  hope  to  identify  a 
subsystem  in  an  aircraft  or  a  subsystem  within  a  subsystem  such  that  selective  maintenance  decisions  can 
be  evaluated  using  the  models  developed  in  this  research. 
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Appendix  A:  Spreadsheet  Instructions  for  Evaluating  Model  P 


Total  Enumeration  Strategy 

File  name:  Total  EnumerationJModelP 

Setting  up  the  spreadsheet  for  use 

When  you  open  the  file,  you  may  receive  a  message  similar  to  the  one  shown  in  Figure  A  1. 


^  . ^  ..piP"U'i3 

C :  \Dc<umen*.s  arid  t  op\fleet\Modd\ T  old 

certains  macros,  "" 


: ... .. 

' .  •  v  •  ifiA 


Macros  may  contab  viruses.  It H  always  to  doable  macros,  but  If  the  :  " 

macros  are  legitimate,  you  might  kt$e  some  functionary; 


if  Enable  Macros  j  Jjorelnfo 


Figure  A  1 :  Macro  Notification 

All  the  Visual  Basic®  code  used  to  evaluate  the  model  is  written  within  a  macro.  Therefore,  you  should 
click  Enable  Macros.  This  will  open  the  Inputs  worksheet  shown  in  Figure  A  2.  Next,  fill  in  the  number 
of  systems,  q,  and  subsystems,  m,  in  the  scenario  you  wish  to  evaluate.  Then,  determine  the  number  of 
limited  maintenance  resources,  5.  (Note  that  for  problems  involving  more  than  4  systems,  computation 
times  become  extremely  long.  Therefore,  large  scenarios  should  be  evaluated  using  the  GA.)  Once  you 
have  entered  these  values,  activate  a  blank  cell  (by  clicking  it),  and  click  the  “Reset  Fields”  button.  Figure 
A  3  shows  the  updated  Inputs  worksheet  for  the  scenario  in  which  q  =  2,m=z3,  and  5  =  2. 
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Figure  A  2:  Inputs  Worksheet 
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Figure  A  3:  Updated  Inputs  Worksheet 


Entering  model  parameter  values 

For  each  scenario  you  wish  to  evaluate,  you  must  enter  the  model  parameters.  The  first  parameters  of 
interest  are  the  reliabilities  of  components  in  each  subsystem,  rh  and  the  number  of  independent  and 
identical  copies  of  components  in  each  subsystem,  rij.  These  values  are  entered  as  shown  in  Figure  A  4. 
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Figure  A  4:  Input  Parameters  for  Component  Reliability  and  Number  of  Components 


Next,  you  must  enter  the  amount  of  resource  consumed  by  repairing  a  type  j  component,  «//  and  the 
amount  of  each  resource  available,  /?,.  These  values  must  be  entered  for  each  of  the  s  limited  maintenance 
resources  (in  this  case,  s  =  2)  as  shown  in  Figure  A  5. 
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Figure  A  5:  Input  Parameters  for  Limited  Maintenance  Resources 

Consider  the  numerical  example  from  Section  3.  Recall  that  q  =  2,  m  =  3,  and  5  =  2,  and  note  that  (:l\  =  1 6 
and  f3i=  10.  The  remaining  model  parameters  are  shown  in  Table  A  1.  Figure  A  6  provides  a  snapshot  of 
the  Inputs  worksheet  with  the  appropriate  parameter  inputs. 


Table  A  1:  Model  Parameters  for  Example  Set  of  Systems 
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Figure  A  6:  Inputs  Worksheet  with  Example  Parameter  Values 


Running  an  experiment 

Once  you  have  entered  all  the  model  parameters,  you  are  now  ready  to  run  an  experiment.  Determine  the 
number  of  failed  components  in  each  subsystem.  If  we  again  consider  the  numerical  example  from 
Section  3,  the  failed  components  for  this  scenario  are  shown  in  Table  A  2.  These  values  are  entered  into 
the  spreadsheet  in  the  same  manner  as  the  other  model  parameters. 
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To  run  the  model,  simply  activate  a  blank  cell  and  click  New  Go.  You  will  be  shown  the  total  number  of 
feasible  solutions  evaluated,  the  run  time,  and  the  optimal  solution  in  Figure  A  7.  For  this  scenario,  72 
solutions  were  evaluated  in  less  than  one  hundredth  of  a  second.  The  optimal  solution  is  to  repair  one 
component  in  each  subsystem  of  system  1  (i.e.  du  =  d]2  =  d\$  =  1)  and  one  component  in  subsystems  2 
and  3  of  system  2  (i.e.  d2\  =  0,  d22  =  d2i  =  1).  The  resulting  overall  reliability  is  0.9479.  The  solution  of  the 
model  is  also  shown  on  the  Output  worksheet  as  shown  in  Figure  A  8. 
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Figure  A  7:  Model  Output 
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Figure  A  8:  Output  Worksheet 
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Genetic  Algorithm 
File  name:  GA_ModelP 
Setting  up  the  spreadsheet 

When  you  open  the  file,  you  may  receive  a  message  similar  to  the  one  shown  in  Figure  A  9. 


Figure  A  9:  Macro  Notification 

All  the  Visual  Basic®  code  used  to  evaluate  the  model  is  written  within  a  macro.  Therefore,  you  should 
click  Enable  Macros.  This  will  open  the  GA  Settings  worksheet  with  the  default  settings.  Recall  that 
population  size  is  the  number  of  candidate  solutions  or  chromosomes  that  the  GA  monitors,  and  the 
number  of  generations  is  the  number  of  times  the  population  undergoes  crossover  and  mutation.  Finally, 
the  seed  value  relates  to  the  random  number  generator.  This  cell  may  have  a  value  from  1  to  100.  This 
worksheet  is  shown  in  Figure  A  10. 


Next,  you  must  enter  the  parameters  of  the  scenario  you  wish  to  evaluate.  This  is  accomplished  by 
entering  data  into  the  Inputs  Worksheet.  Fill  in  the  number  of  systems,  q ,  and  subsystems,  m ,  in  the 
scenario  you  wish  to  evaluate.  Then,  determine  the  number  of  limited  maintenance  resources,  s.  Once  you 
have  entered  these  values,  activate  a  blank  cell  (by  clicking  it),  and  click  Reset  Fields.  Figure  A  11  shows 
the  Inputs  worksheet  for  the  scenario  in  which  q  =  2,  m  =  3,  and  s  -  2. 
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Figure  A  10:  GA  Settings  Worksheet 
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Entering  model  parameter  values 

For  each  scenario  you  wish  to  evaluate,  you  must  enter  the  model  parameters.  The  parameters  of  interest 
are  the  reliabilities  of  components  in  each  subsystem,  /j,  the  number  of  independent  and  identical  copies 
of  components  in  each  subsystem,  w;,  the  amount  of  resource  consumed  by  repairing  a  typey  component, 
dji,  and  the  amount  of  each  resource  available,  /?/.  These  values  must  be  entered  for  each  of  the  s  limited 
maintenance  resources  (in  this  case,  5  =  2)  as  shown  in  Figure  A  12. 

Consider  the  numerical  example  from  Section  3.  Recall  that  q  =  2,  m  =  3,  and  s  =  2,  and  note  that  f3\  -  16 
and  ph.=  10.  The  remaining  model  parameters  are  shown  in  Table  A  3.  Figure  A  12  provides  a  snapshot  of 
the  Inputs  worksheet  with  the  appropriate  parameter  inputs. 

Table  A  3:  Model  Parameters  for  Example  Set  of  Systems 
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Figure  A  12:  Inputs  Worksheet  with  Example  Parameters 


Running  an  experiment 

Once  you  have  entered  all  the  model  parameters,  you  are  now  ready  to  run  an  experiment.  Determine  the 
number  of  failed  components  in  each  subsystem.  If  we  again  consider  the  numerical  example  from 
Section  3,  the  failed  components  for  this  scenario  are  shown  in  Table  A  4.  These  values  are  entered  into 
the  spreadsheet  in  the  same  manner  as  the  other  model  parameters. 

Table  A  4:  Failed  Components  for  Example  Set  of  Systems 


/  £  ■  ■ 

System  (q) 

Subsystem  fQ 

ri  •: 

2  '  it 

3 

1 

an  =  1 

8l2  “  2 

ai3  =  1 
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To  run  the  model,  simply  activate  a  blank  cell  and  click  New  Go.  You  will  be  shown  the  total  number  of 
feasible  solutions  evaluated,  the  run  time,  and  the  optimal  solution,  as  in  Figure  A  13.  For  this  scenario, 
the  GA  ran  for  47  seconds  and  provided  the  same  solution  as  the  total  enumeration  approach.  Therefore, 
the  optimal  solution  is  to  repair  one  component  in  each  subsystem  of  system  1  (that  is,  du  ~  <Jn  —  d\ 3  =  1) 
and  one  component  in  subsystems  2  and  3  of  system  2  (specifically,  d2\  —  0,  d22  =  d2 3  =  1).  The  resulting 
overall  reliability  is  0.9479. 
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Figure  A  13:  Model  Output 
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Appendix  B:  Spreadsheet  Instructions  for  Evaluating  Model 
SP 


File  name:  GA_ModelSP 

Setting  up  the  spreadsheet 

When  you  open  the  file,  you  may  receive  a  message  similar  to  the  one  shown  in  Figure  B  1. 


Mii 

C:\b0curtient5  and  a! 

Is  cootafns  macros. 


MM 


Macros  may  conteh  vi? uses,  It  H  akfsy*  sate  to  disabte  maoos>  but  if  the 
macros  ara  legitimate';  you  might  k&e  som  functionary.  . 


•s  Macros 


vgpnr  Info  :: 


Figure  B  1 :  Macro  Notification 

All  the  Visual  Basic®  code  used  to  evaluate  the  model  is  written  within  a  macro.  Therefore,  you  should 
click  Enable  Macros.  This  will  open  the  GA  Settings  worksheet  with  the  default  settings.  Recall  that 
population  size  is  the  number  of  candidate  solutions  or  chromosomes  the  GA  monitors,  and  the  number  of 
generations  is  the  number  of  times  the  population  undergoes  crossover  and  mutation.  Finally,  the  seed 
value  relates  to  the  random  number  generator.  This  cell  may  have  a  value  from  1  to  100.  This  worksheet 
is  shown  in  Figure  B  2. 


■11 


ssm 


fgrflvst  Jodi  « qpWton t»  twfr  *  w  0  X 

Ar*  -10  *iB  /  U  M  *  SI  $  %  lf~l  Z  *  &  *  & 

vO  S  B  ^  '  ’i  R  V  ?  &  *  <f  '  r  *  $*  V,»  it  R  100%  -  05  „ 


f-ta  $%  $r i4i^a 


B? 


. b  ~i  r 


0 


G  i  H 


J  j . ' 

l  ?  j  Popiilutlon  Size  32 
'  J  j  Generation*  1003 ... 

■  4  J  Seed  Value  1 

$  | 

J4  . . . j . . . •  . . . : 

hi  V  SeM ^ 

•:::p:igjsdy  L,;,4Tv:;.:  •  " 


.  ,j  '  »ir 

«.W  _ d 


Figure  B  2:  GA  Settings  Worksheet 
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Next,  you  must  enter  the  parameters  of  the  scenario  you  wish  to  evaluate.  This  is  accomplished  by 
entering  data  into  the  Inputs  Worksheet.  Fill  in  the  number  of  systems,  q,  and  subsystems,  m,  in  the 
scenario  you  wish  to  evaluate.  Then,  determine  the  number  of  limited  maintenance  resources,  s.  Once  you 
have  entered  these  values,  activate  a  blank  cell  (by  clicking  it),  and  click  Reset  Fields. 

Entering  model  parameter  values 

For  each  scenario  you  wish  to  evaluate,  you  must  enter  the  model  parameters.  The  parameters  of  interest 
are  the  failure  rates  of  components  in  each  subsystem,  Ap  the  number  of  independent  and  identical  copies 
of  components  in  each  subsystem,  rip  the  amount  of  resource  consumed  by  repairing  a  type  j  component, 
a,i,  and  the  amount  of  each  resource  available,  /?/.  These  values  must  be  entered  for  each  of  the  s  limited 
maintenance  resources  (in  this  case,  s  =  2)  as  shown  in  Figure  B  3. 

Consider  the  numerical  example  from  Section  4.  Recall  that  q  =  2,m  =  3,  and  5  =  2,  and  note  that  j3\  =  16 
and  =  14.  The  next  mission  durations  (in  hours)  are  t\  =  4,  t2  =  1.75,  tj,  =  2,  U  ~  0.5,  and  the  remaining 
model  parameters  are  shown  in  Table  B  1.  Figure  B  3  provides  a  snapshot  of  the  Inputs  worksheet  w'ith 
the  appropriate  parameter  inputs. 


Table  B  1:  Model  Parameters  for  Example  Set  of  Systems 


Subsystem 

Aj  (in  failures  per  hour) 

ajr 

Ofl 

1 

B 

0.10 

2 

3 

2 

a 

0.16 

3 

1 

1  3 

2 

0.06 

2 

2 
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Figure  B  3:  Inputs  Worksheet  with  Example  Parameters 


Running  an  experiment 

Once  you  have  input  all  of  the  model  parameters,  you  are  now  ready  to  run  an  experiment.  Determine  the 
number  of  failed  components  in  each  subsystem.  If  we  again  consider  the  numerical  example  from 
Section  4,  the  failed  components  for  this  scenario  are  shown  in  Table  B  2.  These  values  are  entered  into 
the  spreadsheet  in  the  same  manner  as  the  other  model  parameters. 
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Table  B  2:  Failed  Components  for  Example  Set  of  Systems 


System  (g) 

Subsystem  (i)  %> 

2 

yl'Z 
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an  =  1 

a12  =  1 

9i3  =  1 
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^21  =  1 

322  =  2 

a23  —  0 

3 

^31  ~  1 

a32  =  2 

933  “  1 

4 

941  “  0 

942  ~  0 

^43  =  1 

To  run  the  model,  simply  activate  a  blank  cell  and  click  New  Go.  For  this  scenario,  the  GA  ran  for  one 
minute  and  eight  seconds,  and  the  output  shown  on  the  first  line  of  the  Print  Outs  worksheet.  As  you  can 
see  in  Figure  B  1,  the  GA  identified  two  different  solutions  having  an  overall  reliability  of  0.5347. 
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Figure  B  4:  Model  Output 
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Appendix  C:  Spreadsheet  Instructions  for  Evaluating  Model 
DSP 

File  name:  GA_ModelDSP 

Setting  up  the  spreadsheet 

When  you  open  the  file,  you  may  receive  a  message  similar  to  the  one  shown  in  Figure  C 1. 


Microsoft  iKOlt 


. 

CAProumente  mi  Sett ©1 
Enum^eta  macros. 


Haoor-  may  contain  viruses.  It  &  always  safe  to  disable  macrosjiiH-f  ;the 
macros  are  fegPjm-atts,  you  b mmm-  furtttbnafi^* 


JgfeaSe J  InabbMaaos  More  Info' 


Figure  C  1:  Macro  Notification 

All  the  Visual  Basic®  code  used  to  evaluate  the  model  is  written  within  a  macro.  Therefore,  you  should 
click  Enable  Macros.  This  will  open  the  GA  Settings  worksheet  with  the  default  settings.  Recall  that 
population  size  is  the  number  of  candidate  solutions  or  chromosomes  that  the  GA  monitors,  and  the 
number  of  generations  is  the  number  of  times  the  population  undergoes  crossover  and  mutation.  Finally, 
the  seed  value  relates  to  the  random  number  generator.  This  cell  may  have  a  value  from  1  to  100.  This 
worksheet  is  shown  in  Figure  C  2. 
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Figure  C  2:  GA  Settings  Worksheet 


Next,  you  must  enter  the  parameters  of  the  scenario  you  wish  to  evaluate.  This  is  accomplished  by 
entering  data  into  the  Inputs  Worksheet.  Fill  in  the  number  of  systems,  q ,  and  subsystems,  m ,  in  the 
scenario  you  wish  to  evaluate.  Then,  determine  the  number  of  limited  maintenance  resources,  s.  Once  you 
have  entered  these  values,  activate  a  blank  cell  (by  clicking  it),  and  then  click  Reset  Fields. 

Entering  model  parameter  values 

For  each  scenario  you  wish  to  evaluate,  you  must  enter  the  model  parameters.  The  parameters  of  interest 
are  the  failure  rates  of  components  in  each  subsystem,  Xh  the  number  of  independent  and  identical  copies 
of  components  in  each  subsystem,  nh  the  amount  of  resource  consumed  by  repairing  a  type  j  component, 
a//,  and  the  amount  of  each  resource  available  in  each  time  bucket,  /?/.  These  values  must  be  entered  for 
each  of  the  s  limited  maintenance  resources  (in  this  case,  s  =  2)  as  shown  in  Figure  C  3 . 

Consider  the  numerical  example  from  Section  5 .  Recall  that  q  =  2,  m  =  3,  and  5  =  2,  and  note  that  f3\  =  8 
and  =  6.  The  next  mission  durations  (in  hours)  are  t\  =  4,  h  =  1 .75,  1 3  =  2 ,  /4  =  0.5,  and  the  remaining 
model  parameters  are  shown  in  Table  C  L  Figure  C  3  provides  a  snapshot  of  the  Inputs  worksheet  with 
the  appropriate  parameter  inputs. 
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Table  C  1:  Model  Parameters  for  Example  Set  of  Systems 


Subsystem 

nJ 

Aj  (in  failures  per  hour) 
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Figure  C  3:  Inputs  Worksheet  with  Example  Parameters 
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Running  an  experiment 

Once  you  have  entered  all  the  model  parameters,  you  are  now  ready  to  run  an  experiment.  Determine  the 
number  of  failed  components  in  each  subsystem.  If  we  again  consider  the  numerical  example  from 
Section  5,  the  failed  components  as  well  as  the  return  times  for  current  mission  and  start  times  for  future 
missions  are  shown  in  Figure  C  4.  These  values  are  entered  into  the  spreadsheet  as  shown  in  Figure  C  5. 


Figure  C  4:  Parameters  for  Dynamic  Scenario  Example 
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To  run  the  model,  simply  activate  a  blank  cell  and  click  New  Go.  The  solution  for  the  GA  is  shown  in  the 
Output  worksheet.  For  this  scenario,  the  GA  ran  for  195  seconds.  The  suggested  maintenance  actions  and 
mission  assignments  are  shown  in  Figure  C  6. 
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Figure  C  6:  Suggested  Maintenance  Actions  and  Mission  Assignment 
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